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Abstract. We investigate the star-formation history of gas-rich dwarf galaxies, 
taking account of the dynamical evolution of their neutral hydrogen (Hi) enve- 
lope. Gas-rich dwarfs are classified into blue compact dwarfs (BCDs) and dwarf 
irregulars (dlrrs). In this paper, their Hi envelope is clearly shown not to be 
blown away by their stellar feedback. This is concluded since the observed star- 
formation rate (SFR) of gas-rich dwarfs is generally smaller than a critical SFR, 
f/jcrit, at which stellar feedback accelerates the Hi envelope to the escape velocity. 
From this standpoint and the chemical property of sample BCDs, we suggest two 
possibilities; (1) The Hi gas in the envelope of BCDs is consumed to fuel their 
star-formation; and (2) BCDs have a similar star-formation history. We also dis- 
cuss morphological evolution among dwarf galaxies. As long as gas-rich dwarfs 
are isolated, it is difficult for them to evolve into dwarf ellipticals (dEs). When 
the Hi envelope in gas-rich dwarfs is consumed in subsequent star-formation, a 
morphological exchange between BCDs and dlrrs is still expected, consistent with 
previous studies. If the SFR of gas-rich dwarfs was much higher than i/'crit in 
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the past, interestingly, an evolutionary scenario from dEs to gas-rich dwarfs is 
possible. 

Key words. ISM: evolution — ISM: structure — stars: formation — galaxies: dwarf 
1. Introduction 

For a hierarchical cosmic structure formation scenario, the star-formation history of dwarf 
galaxies, which are building blocks of the structure, needs to be known. To determine 
their star- formation history, the effect of stellar feedback (e.g., supernova-driven winds 
and stellar winds) on the dynamical condition of the interstellar medium (ISM) of dwarf 
galaxies is essential, since their gravitational potential well is rather shallow (Skillman 
& Bender 1995 for a review). Thus, the dynamics of the ISM in dwarfs is a central issue 
for their evolution. In this paper, we study isolated dwarf galaxies to investigate their 
evolution, and examine the effect of stellar feedback on the star-formation rate (SFR) 
(Larson 1974; Saito 1979; Dekel & Silk 1986). 

Dwarf galaxies whose star-formation is currently observable have sufficient amounts 
of Hi gas, compared to their dynamical mass (e.g. Thuan & Martin 1981). These gas- 
rich dwarfs are categorized as dwarf irregular galaxies (dlrrs) and blue compact dwarf 
galaxies (BCDs) by their SFR relative to their size. Since their Hi gas mostly surrounds 
their central star- forming regions (e.g. van Zee et al. 1998), we call it an Hi envelope in 
this paper. Although it is widely believed that the Hi envelope is gravitationally bound to 
host galaxies, the dynamical evolution of the Hi envelope is not completely understood. 
This is not only because the peaks of the rotation curves in the Hi envelopes are not 
always satisfactorily determined (e.g. Pustilnik et al. 2001), but also because there is 
no confirmation that the Hi gas in the envelope is used for a subsequent star-formation 
(e.g. Thuan 1985). Therefore, it is useful to investigate the dynamical evolution of the 
Hi envelope from another point of view, that is, in the framework of a star formation 
history (see also Legrand et al. 2001). 
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Indeed, it has been suggested that the Hi envelope is related to star-formation activity 
(e.g. Thuan 1985), which results in the morphological exchange among dwarf galaxies 
(e.g. Kormendy 1985). Concerning the morphological evolution, Thuan (1985) suggested 
from photometric observations that gas-rich galaxies like BCDs and dlrrs evolve into 
gas- less galaxies like dwarf ellipticals (dEs). The inverse evolution is also theoretically 
expected: dEs acquire Hi gas from the ambient intergalactic medium (IGM), and evolve 



into dlrrs and BCDs (Silk et al. 1987). According to paito et al. (2000)| , the morphological 
exchange between BCDs and dlrrs is also expected, if Hi gas in the envelope is channeled 
into star- forming regions. Contrary to the above suggestion, it is also shown from the 
difference of their surface brightness (Bothun et al. 1986) and their dynamical (van Zee 
et al. 2001) properties that gas-rich dwarfs do not easily evolve into dEs. Thus, it is still 
premature to judge which model is correct. In this paper, wc investigate this unclear 
star- formation history of gas-rich dwarf galaxies, considering the dynamical condition of 
the Hi envelope which is affected by stellar feedback. 

This paper is organized as follows: In the next section, we construct a dynamical 
model of the Hi envelope from a momentum conservation law. The observational data are 
compared with the model in section 3. The implications for the morphological evolution 
of dwarfs are given in section 4. In the last section, we summarize our conclusions and 
comment on Legrand et al. (2001), whose conclusions are very similar to and consistent 
with those in this paper. 



2. Model description 

In this paper, we try to connect the dynamical condition of the Hi envelope and star- 
formation history of gas-rich dwarf galaxies. For such a complicated problem, a simple 
analytical model is very useful. Indeed, this would allow us to predict the critical limit 
of a SFR where gas-rich dwarfs can sustain the Hi envelope. If a SFR is larger than 
the critical one, the Hi envelope will be blown away by stellar feedback (e.g. Saito 1979; 
Dckcl & Silk 1986). Ferrara & Tolstoy (2000) analyze similar dynamical conditions of 
Hi envelopes in dwarf galaxies with their dynamical mass being blown-away (all gas 
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is removed), blown-out (partial gas is removed), or no-mass- loss (no gas is removed) 
(classification proposed by De Young & Heckman 1994). These hydrodynamical processes 
are numerically examined in Mac Low & Ferrara (1999), and their results are confirmed 
by Silich & Tenorio-Tagle (2001). Silich & Tcnorio-Tagle (2001) have also shown that 
the low density halo suppresses the outflow of gas from the main body of a galaxy. From 
this standpoint, the circulation of the ISM is expected to be like a chimney (Kunth & 
Ostlin(2000) for a review). This means that a closed box approximation is adequate to 
consider the evolution of gas-rich dwarfs whose ISM is not blown out or blown away. To 
confirm these results with a SFR that can be easily observed compared to dynamical 
mass (Pustilnik et al. 2001), we reformulate a critical value that determines a condition 
for the Hi envelope. Consequently, as shown in the following, we can get intuitive pictures 
of the star-formation history and dynamical properties of the Hi envelope. We have newly 
considered this useful representation which is similar to that of Legrand et al. (2001). 

Before discussing the critical SFR, we review the stellar feedback (i.e., supernovae 
and stellar winds) effect on the Hi envelope. Stellar feedback causes turbulence in the 
ISM, and can make super-bubbles. Indeed, turbulence in the Hi envelope is observed 
(van Zee et al. 1998), and super-bubbles are also observed in some BCDs (e.g. Gil de 
Paz et al. 1999). In a BCD, the typical size of a super-bubble is larger than that of 
its star-forming regions, but is smaller than that of its Hi envelope. Martin (1996) thus 
suggested that the Hi envelope must suffer from the stellar feedback. From the fact that 
the Hi envelope is not blown away, we are able to deduce the allowed strength of the 
stellar feedback, the critical SFR, ip c i It (e.g. Silich & Tenorio-Tagle 1998). Legrand et al. 
(2001) have considered a similar condition and developed a self-consistent model that 
treats together both the Hi envelope and the low density halo (including dark matter). 
On this point, our discussion can provide further insight into the dynamical evolution of 
the Hi envelope. 

To determine V'crit, we deduce the following dynamical condition between the mo- 
mentum of the Hi envelope, p gas , and that supplied from the total stellar feedback by 
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supernovae, psn (Lamers & Cassinelli 1999), as being 

(1 + A ha lo)Pgas > psN, (1) 

where Ah a \ denotes the contribution from the low density halo. It can reach at most 1 
when a closed box model is satisfied (Legrand et al. 2001), which is a good assumption, 
as confirmed in the next section. This inequality is called "the momentum condition" 
(e.g. Ferrara & Tolstoy 2000). Here, we assume that the feedback energy from stellar 
winds is comparable to that of supernova explosions. 

The momentum of the Hi envelope is expressed as p gas = M gas v gas , where M gas is 
the total mass of the Hi envelope, and v gas is the out-going velocity of the Hi gas. For 
simplicity, we adopt the same velocity for the low density halo to use Eq.(l), which 
does not alter our main conclusion since the escape velocity of the low density halo 
is generally larger than v gas . In order to estimate psn, we consider the evolution of a 
supernova remnant in the Sedov phase. As is discussed in Chevalier (1976), the Sedov 
phase of expanding supernova remnants starts when the swept-up mass of a shell of ISM 
or circumstellar medium is comparable to the ejected mass of a supernova (e.g. Jones, 
Smith, & Straka 1981). Thus, the total momentum supplied from supernova remnants 
during a single generation of star- formation is expressed as psn = V^sn Vg N /v s , where 
Msn is the total mass of supernova remnants during the generation, «sn is the typical 
velocity of the initial expansion of the remnants, and v s is the final expanding velocity 
of the shell. Here, r\ is a coefficient whose value is adopted as 15/77 from Mac Low & 
McCray (1988) and which denotes net energy reduction owing to the shell formation 
and the energy conversion to the thermal energy of the inner coronal gas. Therefore the 
momentum condition is re-written as being 

(l + A halo )M gasVgas > i^^dk. (2) 

We can relate the right hand side of Eq. (2) to a current SFR. Assuming a mass 

function, <p(m), we estimate Msn as 

/ioom 
m(f)(m) dm 

M SN = TToowf ' ( 3 ) 



/J-UUJvj 
m(f)(m) dm 
o.im q 
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where ip is a current SFR in the unit of M© yr _1 , and tsF is the duration of a single 

generation of star-formation in the unit of yr. 

The gas mass of the Hi envelope is known to be dominated by the Hi compared to H2 
(Baronc ct al.2000). This also means that the direct contribution of halo gas to the Hi 
envelope is small. However, we take into account halo gas using a coefficient of (1 + Ah a \ ) 
in Eq. (2), because we consider the whole gas system under a closed box approximation. 
The out-going velocity v gas is assumed to be about the escape velocity of the dwarfs with 
a fiat rotation curve, given as 

/ 2G M gas , /Mgas(l + A^)\ 2/7 ,-2,7(F\- 2/ \ , 

where we adopt the Burkert profile of gas-rich dwarf galaxies (Burkert 1995) using a 
radius, i? ga i = 6.9 x 10 _4 M gas 3 / 7 F~ 3 / 7 kp C . Here, G is the gravitational constant, and 
F is the gas mass ratio to dynamical mass of a galaxy. The typical parameters are taken 
from the data of BCDs in Table 1. It should be noted that the value of F used is a 
larger than F in Burkert (1995), 0.2, and thus that our underestimated Vcrit reinforces 
all gas-rich dwarfs. 

Equating both sides of the momentum condition, we obtain ^crit once a mass function 

has been determined. The Salpeter mass function, <p(m); <p(m) oc to -2 ' 35 (0.1 < to < 

100 M Q ) (Salpeter 1955) is assumed here. Hence, 

9/7 

(l + A halo )- 2/7 M yr-\ (5) 

where we adopt the typical values t$F — 10 7 yrs, F = 0.5, A'lg as = 10 8 M , vsn = 1000 
km s _1 , and v s = 100 km s _1 . A correction coefficient of C cs f is for the component of 
continuous star- formation, C cs f = £ a geV'contAsFV , crit, while V'cont is a continuous SFR 
and t age is galaxy age. Here C cs f is generally less than 1/10 the size of the component of 
current star-formation (Legrand et al. 2001; their table 2). Thus, we neglect it since we 
are interested in the current condition of the Hi envelope. As long as the observed SFR 
is well below Vcrit, the Hi envelope is not blown away by stellar feedback. If we adopt 
another IMF and another set of mass ranges, there is generally 20 percent uncertainty 
(e.g., Inoue et al. 2000). 



(1 + Ccsf^crit = 1.1 X 



M gas (1 + A ha lo) 

10 8 M Q 



Y.Y.Tajiri and H.Kamaya: Neutral Hydrogen of BCD 7 
In the next section, we compare i/'crit with the SFR of dwarfs that is estimated from ob- 
servations. For this purpose, we deduce f/> C rit with luminosity of Ha, Ljja (e.g. Kennicutt 
1983), using the relation ijj = 2.32 x 10~ 42 Ln a M & yr _1 that is compiled in and used by 
Legrand et al. (2001) and based on Leitherer & Heckman (1995). Thus, we can obtain a 
critical Lhq from eq.(^J): 

£Ha,crit = 4.1 x 10 



9/7 

(l + A halo ) 9/7 erg. (6) 



10 8 M, 

3. Observational investigation 

3.1. Hi envelope and star-formation rate 

We compare the observational data of Lh q with eq. (||) . The sample galaxies are selected 
from |van Zee ct al. (1998)| , |Sage et al. (1992)| , |van Zee (2000] , and [van Zee (20"01)| , ac- 



cording to the condition that they are isolated systems. Regarding four galaxies found in 
the several papers, the latest data of each galaxy are adopted. We also plot dlrrs in the 



primary sample of van Zee (2000) and van Zee (2001) to present how well our criterion is 
satisfied by dlrrs as well as BCDs. The sample BCDs are summarized in Table 1. The re- 
sult is depicted in Figure 1 . Filled circles are BCDs and the pluses are dlrrs. We draw four 
model lines. The main line is the thick solid line from eq.@ with Ahaio = 0.0. Clearly, all 
the sample BCDs are located below this line. That is, all the Hi envelopes are not blown 



away by the current stellar feedback. This result is the same as Ferrara fc Tolstoy (2000) , 
since the dynamical mass-range of our sample galaxies is from ~ 10 7 M Q to ~ 10 9 Af Q . 
The thick dashed line is from eq. (||) with average internal extinction of gas-rich dwarf 
galaxies (E(B-V)=0.2) (Hunter & Gallagher 1986) and Ahaio = 0.0. Even if wc consider 
the effect of extinction, the conclusion is not altered. The thin lines are the cases with 
^haio = 1-0. The thin lines are also sufficiently above the sample BCDs. This strongly 
indicates that the Hi envelope is not blown out by the support of the low density halo 
as predicted in Silich & Tenorio-Tagle (1998). 

Wc find one interesting sample galaxy, Mrk 67, which is plotted at (M(Hl),L(Ha)) = 
(1O.O 7 4 M , 10.0 39 9 ) in Figure 1. It is a Wolf-Rayet galaxy (Conti 1991). This galaxy is 
located near the solid line. The Hi envelope is strongly affected by stellar feedback. Thus, 
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galaxy 


logM(HI)[M ] 


logL(Ha)[erg] 


F 


12+log(0/H) 


[Fc/O] 


Ref 


IIZw40 


8.64 


41.23 


0.68 


8.13 




a,d 


UGC4483 


7.57 


38.64 


0.32 


7.52 




a 


UM439 


8.49 


40.07 


0.36 


8.05 




a 


UM461 


8.23 


40.00 




7.74 




a 


UM462 


8.42 


40.52 


0.51 


7.98 




a 


Haro2 


8.68 


41.00 


0.31 


8.4 




b,d 


Haro3 


8.76 


41.03 


0.47 


8.3 




b,d 


UM465 


7.71 


39.90 




8.9 




b 


Mrk67 


7.36 


39.95 


0.20 


8.09 




b,d 


Mrk900 


8.18 


40.37 




8.5 




b 


Mrk328 


8.23 


40.51 


0.09 


8.5 




b,d 


UGC11755 


8.20 


39.99 








e 


UGCA439 


8.52 


39.87 


0.66 






d,e 


Mrk600 


8.52 




1.51 


7.83 


0.17 


c,d 


Mrk5 


8.11 




0.50 


8.04 


0.14 


c,d 


Mrk71 


9.00 




0.68 


7.85 


0.49 


c,d 


IZW18 


7.81 




2.09 


7.18 


0.38 


c,d 


Mrk22 


8.18 




0.69 


8.00 


0.57 


c,d 


Mrk36 


7.28 




0.79 


7.81 


0.37 


c,d 


VIIZw403 


7.58 




0.61 


7.69 




c,d 


Mrk750 


7.23 




0.18 


8.11 


0.50 


c,d 


Mrk209 


7.70 




1.13 


7.77 


0.49 


c,d 


Mrk59 


8.66 




0.36 


7.99 


0.50 


c,d 



Table 1. BCD sample galaxies. The date are taken from a (van Zee et al. 1998), b (Sage 
et al. 1992), c (Izotov & Thuan 1999), d (Thuan & Martin 1981) and e (Van Zee 2000) 
and (van Zee 2001). Blanks mean no data. 

if we want to know the precise effect of the low density halo, the detailed observation of 
the more extended Hi envelope of Mrk 67 is important. 
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3.2. Star-formation history of BCD 

In this paper, we are interested in the star-formation history of isolated BCDs, which 
can be given by their chemical abundance in their star-forming region. We present this 
in Figure 2, while the data of sample galaxies is listed in Table 1. In panel (a), we plot 
the metallicity of BCDs and the gas mass to the dynamical mass ratio. The metallicity 
is determined in the current star-forming region (van Zee et al. 1998; Izotov & Thuan 
1999). The figure shows that the mass ratio becomes smaller as the metallicity grows 
larger. Thus, panel (a) supports a closed box assumption for the isolated BCDs, as is 
originally found in Lcqueux et al. (1979). For a detailed review, see Kunth & Ostlin 
(2000). Furthermore, as long as a closed box assumption is satisfied, it can suggest that 
the Hi gas in the envelope is used for the subsequent star-formation in the current star- 
forming region. This is consistent with Figure 1, as long as Hi gas is generally bound to 
the host galaxies. It should be noted that higher mass ratios than 1 are obtained simply 
from observational uncertainties of Hi gas and/or estimation of dynamical mass (Thuan 
& Martin 1981). Panel (b) represents the metal abundance and the same mass ratio. 
This is also interesting since there are no significant differences in abundance among the 
BCDs (The abundance is located around 0.4, which is a typical value of SN n). Panel 
(b) suggests that the star-formation histories of BCDs are not very different from each 
other, as discussed in Kinman & Davidson (1981). Thus, we can conclude that almost 
all BCDs have experienced a similar star-formation history. 

4. Discussion 

Our results can provide some insight into a morphological evolutionary model among 
dwarf galaxies (Kormendy 1985). The most standard model for this kind of research, 
which discusses the Hi envelope, is Thuan (1985). They predict that gas-rich dwarfs 
evolve into dEs after the Hi envelope is stripped for some reasons, such as: (1) interaction 
with clusters and/or galaxies; and (2) stellar feedback. Since the gas-rich dwarf sample in 
this paper is isolated, only the stellar feedback effect is expected to be important for the 
morphological exchange. Importantly, the current stellar feedback does not blow away 
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the Hi envelope as shown in Figure 1. Presumably, it does not blow out the Hi envelope 
because of the effect of the low density halo (see also Legrand et al. 2001). Thus, we 
conclude that the isolated gas-rich dwarfs can keep their morphological type, until they 
no longer experience any interactions. As long as the the Hi envelope has been retained, 
the kinematics of Hi envelope is expected to be different from that of stars in dEs. This 
is checked by estimating the specific angular momentum, and indeed a difference exists 
(van Zee et al. 2001). Thus, our conclusion is consistent with other studies, since isolated 
BCDs never evolve into dEs. 

In Figure 2, we suggest that the Hi gas is used for the subsequent star- formation. This 
might permit BCDs and dlrrs to evolve to dEs. However, the consumption time-scale, 
which is defined as the mass of the Hi envelope divided by the current star-formation rate, 
should be longer than 10 10 years. It is difficult for current BCDs and dlrrs to evolve into 
dEs. Our discussion may confirm the observational implication on its difficulty (Ferguson 
& Binggcli 1994) in terms of a SFR. We comment that the photoionization effect can be 
important during the evolution from dlrrs to dEs (Ferrara & Tolstoy 2000), while this 
effect cannot be examined in our framework. This possibility, related to photoionization 
effects, will be discussed in a future work. 

There is another possibility for the morphological change between dlrrs and BCDs. 
According to Saito et al. (2000), gas-rich dwarfs change their morphology between BCDs 
and dlrrs. This means that the star-formation activity is related to the extent of the 
Hi envelope. Here, we re-consider this evolution in Figure 1: dlrrs with massive Hi gas 
move upper-leftward in Figure 1 and are recognized to be BCDs, when the Hi gas falls 
onto galaxies and is consumed for the successive star-formation period. On the contrary, 
when BCDs decrease their star-formation, they move downward in Figure 1, and are 
recognized as dlrrs. Thus, our plot of Figure 1 is consistent with the preliminary idea of 
Saito et al.(2000). The trend of consumption of Hi mass in Figure 2(a) may support this 
evolutionary model of gas- rich dwarfs. 

An inverse evolutionary scenario from dEs to gas-rich dwarfs has been suggested by 
Silk et al. (1987). The IGM was heated and enriched by galactic winds at higher redshift. 
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When galaxy groups formed recently, this gas was compressed, cooled and accreted onto 
dwarf galaxies with some metallicity. Thus, dEs evolve into dlrrs and BCDs. This scenario 
is re-interpreted if we exchange the IGM for the low density halo around the BCDs. Since 
the low density halo is expected to be bound to the galaxy, a closed box assumption can 
also be invoked here. Furthermore, both IGM and the low density halo have a similar role 
in the dynamical evolution of Hi envelope (Silich & Tenorio-Tagle 1998). Also, if a past 
SFR was larger than Vcrit, the Hi envelope could have been blown away from the host 
dwarf galaxy, while being retained in the galaxy by the low density halo. Hence, when 
future observations will reveal the large past SFR of BCDs, unlike current observations 
of their spectral energy distribution, then, it will be interesting for us to re-examine the 
possibility of the morphological evolution model of Silk et al. (1987) with caution. 

5. Summary 

In this paper, we have shown that the Hi envelope are not easily blown away by stellar 
feedback without the low density halo but with dark matter, since the envelope is not 
accelerated to the escape velocity (see also Legrand et al. 2001). This dynamical state is 
observationally confirmed for isolated gas-rich dwarfs. Furthermore, we have confirmed 
the possibility that all the isolated BCD samples have a similar star-formation history, as 
is suggested from their metallicity and Hi mass consumption. Our results support that 
(1) It is difficult for the isolated gas-rich dwarfs to evolve into dEs, since the present 
SFRs of gas-rich dwarfs are not so large. They will change their morphological types 
only when they suffer exceptionally intense star-burst and/or interaction among galaxies, 
and/or clusters of galaxies. (2) The morphological exchange among gas-rich dwarfs may 
be possible. But, as a future work, we must describe the star-formation mechanism due to 
the Hi gas fueled into star-forming regions, suggested from the anti-correlation between 
Hi mass and metallicity. The effect of the low density halo may be essential. (3) A 
morphological evolutionary model from dEs to gas-rich dwarfs is possible if the Hi gas falls 
onto the host galaxy. But, even in this case, the low density halo should be investigated 
further. Recently, Legrand et al. (2001) has reported a very similar result to our Figure 
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1. The essential meaning of the two results are the same and consistent with each other, 
such that the ISM is sustained in the BCDs themselves. However, we consider especially 
the Hi envelope, and find that the Hi envelope is neither blown away nor blown out from 
BCDs by means of the current stellar feedback. 
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Fig. 1. The luminosity of Ha and the mass of Hi gas instead of M gas . The filled circles 
are BCDs and the pluses are dlrrs. The solid lines are from eq. (|^) and the dashed lines 
are from eq. (^) with average internal extinction in dwarf galaxies (E(B-V)=0.2). The 
thick and thin lines are from eq. (p) at Ahaio = and Ahaio = 1, respectively. 
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Fig. 2. a) The anti-correlation between the metallicity and the ratio of Hi gas mass to the 
dynamical mass of samples in Table 1. The filled circles denotes samples with the mass 



ratio from Thuan & Martin (1981), and the plus denotes those of van Zee et al. (1998). b) 
The flat distribution of metallicity abundance to Hi gas mass ratio. The metal enrichment 
is determined at the current star-forming region. 



